Universality of the Kondo effect in quantum dots with ferromagnetic leads 
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We investigate quantum dots in clean single- wall carbon nanotubes with ferromagnetic PdNi-leads 
in the Kondo regime. In most odd Coulomb valleys the Kondo resonance exhibits a pronounced 
splitting, which depends on the tunnel coupling to the leads and an external magnetic field B, and 
only weakly on gate voltage. Using numerical renormalization group calculations, we demonstrate 
that all salient features of the data can be understood using a simple model for the magnetic 
properties of the leads. The magnetoconductance at zero bias and low temperature depends in a 
universal way on gfiB(B — B c ) /k^T^, where 7k is the Kondo temperature and B c the external field 
compensating the splitting. 

PACS numbers: 73.23.Hk, 73.63.Fg, 72.15.Qm, 72.25.-b 



The Kondo effect resulting from the exchange interac- 
tion of a single spin with a bath of conduction electrons 
PQ, is one of the archetypical phenomena of many-body 
physics. Its competition with ferromagnetism and possi- 
ble applications in spintronics [2] have raised wide inter- 
est in the past few years. The Kondo effect in quantum 
dots [5J [4] has, in recent experiments, been investigated 
in the presence of ferromagnetic (FM) leads [5HZ]- It 
was found that the Kondo resonance, usually observed 
at zero bias in the odd Coulomb blockade (CB) valleys, 
is split into two peaks at finite bias [5] . The splitting con- 
sists of a term depending logarithmically on gate voltage 
[6j[7], and, as demonstrated here, a second term nearly 
independent of gate voltage. These phenomena were pre- 
dicted theoretically [8HTT], attributing the splitting of 
the Kondo resonance to a tunneling induced exchange 
field, which results from the magnetic polarization of the 
leads. So far no detailed and quantitative comparison of 
the measured conductance with the theory has been un- 
dertaken to verify whether the simplistic description of 
FM leads used in Refs. I8HTT1 has quantitative predictive 
power. The latter would be needed for future spintron- 
ics applications that exploit the lead-induced local spin 
splitting, e.g., spin filtering. 

In this Letter we present low-temperature transport 
measurements of a single wall carbon nanotube quantum 
dot with PdNi leads. We concentrate on the less studied 
gate-independent contribution of the exchange splitting 
of the Kondo resonance and attribute it to the saturation 
magnetization of the contact material. We show that 
the evolution of the conductance with magnetic field and 
gate voltage can be understood within a simple model for 
the magnetization and polarization in the FM leads, by 
presenting numerical renormalization group (NRG) cal- 
culations for this model, using parameters extracted from 
experiment. Moreover, by comparing resonances of dif- 
ferent transparency, we demonstrate a universal scaling 
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FIG. 1. (Color online) Differential conductance dI/dV s d ver- 
sus bias voltage V^d and gate voltage lg a te at T — 25 mK and 
B — T. The CB regions are numbered for future reference. 

of the magnetic field dependence of the Kondo conduc- 
tance, proving that the magnetization of the leads can 
indeed be viewed as an exchange field, which acts analo- 
gously to an external magnetic field. 

Experimental setup. — The nanotubes are grown by 
chemical vapor deposition on a highly doped silicon sub- 
strate with a 200 nm thermally grown oxide layer at its 
surface [T2] . The contact electrodes with a thickness 
of roughly 45 nm are subsequently structured by elec- 
tron beam lithography and consist of Pdo.3Ni .7, known 
to generate highly transparent contacts [15] . Transport 
measurements were done in a dilution refrigerator with a 
base temperature of 25 mK. The current I is measured 
in a two point geometry with a voltage bias V^d applied 
to the source contact. 

The differential conductance G = d7/dV^d is plotted in 
Fig. [T] in color scale, providing the typical charging dia- 
gram of a quantum dot [14] . Our device exhibits regular 
CB oscillations and a clear fourfold symmetry character- 
istic for carbon nanotubes [15] . Coupling to the contacts 
is strong, leading to broad resonance lines and a variety 
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of higher-order processes. Fig. [T] implies the following pa- 
rameters for the nanotube quantum dot jT5j [16] : a charg- 
ing energy U ~ 5meV, a level separation A ~ 9.5 meV, 
and a subband mismatch of about S ~ 1 meV. The tun- 
nel coupling T between leads and quantum dot can be 
inferred from the line width of the conductance peak. 
Between valleys 1 and 2 in Fig. [I] we obtain a full width 
at half maximum (FWHM) of T = 1.1 meV > k B T. 

A striking feature visible in Fig. [I] are lines of enhanced 
conductance at small, approximately constant bias values 
in every second CB diamond. We attribute these lines to 
a spin- 1/2 Kondo conductance anomaly, split into two 
peaks at small opposite bias values due to the presence 
of FM contacts. The peak distance 2As at the center of 
the diamond can be related to a magnetic field scale B — 
As/gfiB via the Zeeman effect. For valley 1 in Fig. [I] this 
gives approximately 2 T for g = 2. Some resonances with 
very low conductance exhibit no measurable splitting. 

Main features of B- dependence: — Figure |g|a) displays 
detailed measurements of valley 1 from FigT^Jfor differ- 
ent values of external magnetic fields almost parallel to 
the nanotube axis. The main observations are: (i) The 
dominant contribution to the splitting is independent of 
VJr a te- (ii) As the field strength increases, the splitting is 
reduced (observed in all investigated cases) until only a 
single apparent peak remains. This field value is referred 
to as compensation field (here B c ~ 2T), since the dom- 
inant gate-independent part of the splitting is compen- 
sated. At higher fields the peak splits again, (iii) Despite 
(i), we observe nevertheless a slight gate dependence, in 
particular near B c . This is most clearly reflected by the 
fact that the touching point of the two resonances moves 
from the left side of the CB diamond for B < B c (cf. 
B = 1.5 T) to the right side for B > B c (cf. B = 2.5 T). 
(iv) The gate dependence increases in strength very close 
to the edges of the CB diamond. 

Tunneling-induced level shifts. — The presence of the 
splitting, its dependence on magnetic field, and a poten- 
tial gate dependence can be consistently explained by the 
renormalization of the quantum dot level energy due to 
charge fluctuations between the dot and the leads. Since 
the density of states (DOS) in the FM contacts is spin 
dependent, the same is true for the charge fluctuations 
and the corresponding energy renormalization. Hence, 
the renormalization results in a spin splitting of the dot 
level Ae = 5e^ — Se ± — g\iBB. For a single impurity An- 
derson model, the correction Se^ (a = t,i) from second 
order perturbation theory is (see, e.g., Ref. [TO]): 
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Here, = T 9/^bB/2 is the quantum level energy 
for spin a, U the charging energy, and T a (u) the spin- 
dependent tunneling rate. f(cj) is the Fermi function. 
From Eq. ([TJ one sees that the splitting is not only 
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FIG. 2. (Color online) Differential conductance versus 
source-drain and gate voltage, in CB region 1 of Fig.[l] for six 
magnetic field values. We compare here (a) experimental data 
with (b) NRG results for the normalized zero-temperature 
spectral function Q — (0)A a (cj), obtained using the 

model and parameters described in the text. n g = 1/2 — Sd/U 
is the dimensionless gate potential. The dashed line in the 
panel for B = 2 T in (a) indicates the gate-dependent contri- 
bution from the polarization for V — 10% (see text). 

a consequence of properties at the Fermi surface, but of 
the full DOS. The first and second terms in Eq. de- 
scribe electron- or hole-like processes, meaning fluctua- 
tions between the states |l,cr) and |2) or |0), respectively 
(the numeral denotes the charge occupation of the quan- 
tum dot). The spin dependent energy corrections Se a are 
negative, as always in second order perturbation theory. 
Consequently, the spin direction that favors fluctuations 
more strongly will have lower renormalized energy. 

Effect of Magnetization. — First we assume a shift be- 
tween bands of equal and constant DOS for different spin 
directions, p^ = p± = po, described by a constant Stoner 
splitting Ast, see Fig. [3j The tunneling induced splitting 
/\ £ {M) ^ ue £ Q i s directly related to the saturation 
magnetization M. = (n^ — n^)/N a = Ast/(2-Do)- Here 
= A) (A) +£f± Ast/2) is the number of spin-a elec- 
trons, £p the Fermi energy, and N a the number of states 
per atom and spin. Starting from Eq. ([T]) we can de- 
duce the spin orientation of the dot ground state as fol- 
lows: Figs. |3jb) and (c) show the phase space available 
(hatched) for quantum charge fluctuations for a spin up 
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FIG. 3. (Color online) Schematic of the level renormalization 
process, (a) We assume flat bands pa{oj) — po with bandwidth 
Do, shifted with respect to each other by a constant Stoner 
splitting Ast- The filling fraction T determines the Fermi 
energy sf- (b) Charge fluctuations for a spin up electron on 
the dot, involving the empty states of the spin up band and 
the occupied states of the spin down band (hatched areas), 
(c) Analogous situation for a spin down electron on the dot; 
the available phase space (hatched) is larger than in (b). 

or a spin down electron residing on the quantum dot, re- 
spectively. Comparing the total sizes of hatched areas for 
(b) and (c), the phase space is larger for the latter, thus 
favoring spin down. Therefore, Ae^ M ^ is always negative, 
meaning that the energy correction due to a finite Stoner 
splitting is larger for the minority spin, i.e. the ground 
state spin is oriented opposite to the magnetization of the 
leads. This explains why the splitting is always initially 
reduced (never increased) when an external field is ap- 
plied (cf. observation (ii) listed above; the magnetization 
direction follows that of the field in our setup). The size 
of this effect depends on the Stoner splitting, i.e., on M. 
For |e<i|j £d + U «Do, Ast, which is compatible with the 
experiment, we obtain, in extension of [TT] , 
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where T = (1 + Sy/Dq)/2 is the filling fraction of the 
band. This shift is independent of the gate voltage (ex- 
plaining observation (i)) because the position sa of the 
level is very close to £p, while the integration in Eq. ([I]) 
ranges over a large fraction of the d-band. 

SQUID measurements allow us to determine a mag- 
netic moment of p = 0.58 p& /atom for our PdNi alloy 
[17], implying an effective magnetization of M = 0.116. 
Ab-initio calculations of the band structure provide a 
filling fraction of T = 0.853 [18 j. Thus, we estimate 
/\ £ ( M ) ~ — 175peV. In transport spectroscopy this 
would lead to conductance peaks split at zero exter- 
nal field by 2As theo /e = 2\Ae^/e\ ~ 350 pV . Given 
the simplicity of our model, this agrees reasonably well 
with the experimentally determined peak distance of 
2Ae exp /e ~ 550 pV. For a more weakly coupled reso- 
nance the predicted peak distance of 60 pV agrees simi- 
larly with the experimentally found value of 105 pV. 

Effect of polarization. — The case of p^ =^ p^ implying 
nonzero polarization V = (p^ — Pi)/(p^ + p\) at £p, has 
already been discussed in Refs. [6j[TT] and earlier publica- 



tions referenced there. Assuming a fiat band with a spin- 
dependent DOS, e.g. p^ > pi but zero M for simplicity, 
quantum charge fluctuations renormalize the quantum 
dot level depending on its position relative to €p. This 
contribution shows a logarithmic divergence for ea — >• 
and £d + U — » [H [11] , resulting in the up- and downward 
bending of the compensated conductance peak towards 
the corners of the diamond (cf. observation (iv)). 

Numerical results. — The quality of our model is re- 
flected by the close correspondence of Figs.[2|a) and (b). 
Fig. [2^b) presents high-quality NRG results for the spec- 
tral function A(uj) [19 versus the normalized gate voltage 
n g , calculated for a single-lead Anderson model with the 
DOS shown in Fig. 3 but with p^ ^ using full density- 
matrix NRG [20j [21] • Using the measured parameters 
of the quantum dot and modeling the ferromagnetism 
in the leads by taking V = 10% and M. = 0.116, very 
good agreement with experiment is found, except for the 
background current at high V s d, which results probably 
from cotunneling processes involving higher levels not in- 
cluded into the model. In the experiment as well as in 
the numerical data, at B = 2 T the gate-independent 
contribution, cf. Eq. is fully compensated and the 
crossing point of the resonances lies in the center of the 
CB diamond. Here, only the (weak) gate-dependent con- 
tribution from V remains (cf. observation (iii)), indicated 
in Fig. [2] {B = 2 T) by a dashed line. By varying V and 
comparing the shape of the conductance peak at B c be- 
tween experiment and theory, we estimate V ~ 10%. 

Universality. — For a quantum dot coupled to normal 
leads, the normalized zero-bias conductance is a univer- 
sal function (1) of T = T/Tk at zero field, and (2) of 
B = gp^B/k^T^ at zero temperature. (We define the 
Kondo temperature T K via G(T = T K )/G(0) = 1/2 at 
zero field). We find, quite remarkably, that both these 
universal features are recovered also for ferromagnetic 
leads, if B is replaced by the effective field SB = B — B c . 
Regarding (1), Fig. Qd) shows that at the compensation 
field, B = B C1 the temperature dependence of the nor- 
malized conductance, G(T)/G(0), agrees with the often- 
used semi-empirical formula g(T) = [1 + (2 X / S — 1)T 2 ] _S , 
with s ~ 0.22 [5- Although the latter behavior is well- 
established for dots coupled to normal leads, its emer- 
gence here is nontrivial: it demonstrates that despite the 
ferromagnetic environment, local spin symmetry can in- 
deed be fully restored by fine-tuning the field to B c . 

The magnetic field dependence (2) has so far attracted 
much less attention [22 . To explore it, Fig. Qa)-(c) 
shows G(B) at fixed T < Tk for several charge states 
differing in Tk, B C1 and Ae(B = 0). The position of the 
conductance peak roughly follows the Zeeman law, with 
slight deviations in the vicinity of the compensation field 
B c [22]. We find that B c and Tk vary independently for 
different charge states, implying different couplings T. 

According to (2), G(B)/G(B C ) should be a univer- 
sal function of 5B = gp^{B — B c ) /k^T^. The lines 
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FIG. 4. (Color online) (a)-(c) Magnetic field dependence of 
the splitting for different charge states. Vertical arrows denote 
the compensation field B c . (d) Scaling plot of G(T) vs. T/Tk 
(symbols), at B — B c , for the charge state shown in Fig.EJc); 
the solid line gives G(0) times the universal function g(T/Tk) 
discussed in the main text, (e) Scaled zero bias conductance 
G(B)/G(B C ) plotted against the effective normalized field SB, 
at fixed T = 50 mK < T K (100 mK for ■). Lines represent 
NRG calculations for several parameter sets. Inset: G(B) vs. 
(B — B c ) before scaling. 

in Fig. |4je) show this curve, calculated by NRG for 
four different sets of model parameters, yielding a good 
scaling collapse. Symbols show experimental data for 
G(B)/G(B C ) vs. SB, for several different gate voltages, 
with Tk extracted by numerical fitting to the NRG re- 
sults. For three data sets taken on the same charge 
state (circles and triangles), scaling works very well and 
agreement with theory is excellent; for the other two 
sets (squares, diamonds), the quality of scaling is re- 
duced at higher SB by an asymmetric background con- 
tribution to the magnetoconductance. Nevertheless, the 
overall agreement between theory and experiment shows 
that the model correctly captures the universal, sample- 
independent features of G{B) as function of SB. 

Conclusions. — We have performed a quantitative 
comparison of the conductance of quantum dots with FM 
contacts, in the Kondo regime, with model NRG calcula- 
tions. The quantitative agreement between experimental 
and numerical data lends strong support to the scenario 
proposed in Refs. |8] and [TTJ the exchange field induced 
by magnetic contacts causes the local level to be split 



by an amount Ae, which adds a constant offset to the 
Zeeman splitting induced by an external magnetic field. 
When this offset is compensated by a suitably fine-tuned 
field 5 C , universal scaling features of the Kondo effect are 
recovered. The control of Ae(B) through the selection of 
materials with proper M. and V (high or small M. opti- 
mize peak splitting or gate tunability, respectively) may 
prove useful for future applications in spintronics. 
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